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PROCESS FOR USING SURFACE STRAIN 
MEASUREMENTS TO OBTAIN 

OPERATIONAL LOADS FOR COMPLEX 
STRUCTURES 

5 

STATEMENT OF GOVERNMENT INTEREST 

The invention described herein may be manufactured and 
used by or for the Government of the United States of 
America for governmental purposes without payment of any io 
royalties thereon or therefor. 

BACKGROUND OF THE INVENTION 

1 . Field of the Invention 15 

The present invention relates generally to processes to 

obtain externally applied, out-of-plane operational load data 
using surface in-plane strain measurements on complex 
structures, particularly to real-time methods to calculate 
operational loads for complex structures, and most particu- 2Q 
larly an ultra-efficient process of calculating externally 
applied bending and torsional operational loads in real-time 
on complex structures using strain measurements provided 
by fiber Bragg grating technology. 

An important by-product of this process is the efficient 
characterization of out-of-plane bending and torsional stiff- 25 
ness (material properties) of the complex structure. 

2. Description of the Related Art 

Currently, to obtain operational load data for complex 
structures, industry uses a very rigorous and time consuming 
design methodology that relies heavily on computational 30 
methods, such as finite element modeling (FEM). FEM 
requires that the structure to be designed/analyzed be analyti- 
cally discretized into very small linear “elements.” Realistic 
structures are geometrically complex, as are the loads they are 
exposed to and the boundary conditions that restrain them. 35 
Therefore a considerable amount of effort is required to accu- 
rately model the structural response of these realistic struc- 
tures in the loading environment for which they are expected 
to perform. An extraordinary amount of knowledge of the 
structural design details is required a priori (eg. all the local 4Q 
material and geometric properties of the structure, how the 
skins attach to the spars, what the load paths are, the nonlin- 
earity of the joined regions, how the composite design differs 
from the “as built” structure, etc). In addition, the models are 
also only as accurate as the material properties assumed in the 
analysis. These material properties are usually derived from 45 
small (1 inxlO in) uniaxial coupon tests, the results of which 
are assumed to be appropriate for the three-dimensional 
structure. The development of FEM models is extremely 
labor intensive and requires costly experimental validation to 
ensure that the models accurately reflect reality. 50 

More specifically, for real time loads associated with air- 
craft wings, the methodology used today relies upon conduct- 
ing a fairly extensive “strain gage/loads calibration.” This 
method is built upon technology developed in the 1 950s and, 
first, involves the installation of a sparsely distributed number 55 
of conventional strain gage sensors located on the inside of 
the wing structure. An extensive series of tests are then 
required to apply a matrix of independent concentrated loads 
to the wing surface to “calibrate’ the strain gage response to 
each particular load applied, called influence coefficients. 
Software is then used to derive loads equations that can be 60 
used in real time using the strain gage measurements as input. 
The drawback to this method is that it involves a considerable 
amount of labor and cost to perform these strain gage loads 
calibrations from installing strain gage instrumentation, to 
designing and fabricating load system hardware, to assem- 65 
bling the test setup, to conducting the tests, to reducing the 
data and deriving the loads equations. A fairly elaborate load- 
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ing system is required to be designed, fabricated, and 
assembled in order to conduct the tests. 

Due to the shortcomings of the current real time, opera- 
tional loads analysis discussed above, it is desired to provide 
an improved, real-time, method to calculate operational load 
data for complex structures, and, more particularly aircraft 
wings, that significantly reduces time and cost. 

SUMMARY OF THE INVENTION 

The invention is an improved process for using in-plane 
surface strain data to obtain real-time, externally applied, 
out-of-plane operational bending and torsional loads data for 
complex structures that significantly reduces the time and 
cost versus current methods. 

Accordingly, it is an objective of this invention to provide 
real-time, externally applied, out-of-plane operational loads 
data (bending and torsion) for complex structures. 

It is another objective of this invention to reduce the time 
necessary to obtain real-time, externally applied, out-of- 
plane operational loads data (bending, and torsion) for com- 
plex structures. 

It is a further objective of this invention to reduce the cost 
associated with obtaining real-time, operational loads (bend- 
ing and torsion) for complex structures. 

This invention meets these and other objectives related to 
determining the real-time, operational loads on complex 
structures by providing an improved method of obtaining said 
data that significantly reduces the cost and time associated 
with the current analysis. 

In general, the invention comprises installing a very fine, 
spatially -refined grid of fiber optic sensors on the complex 
structure to be analyzed. For bending, a single-point load is 
externally applied to the complex structure at the free end to 
obtain the bending properties of the complex structure to the 
“as built” design. Similarly, a single force-couple (moment) is 
applied at the free end of the structure to obtain torsional 
structural properties of the complex structure. Because the 
local bending and torsional stiffness properties of the com- 
plex structure can be to determined at every fiber optic sensor 
location, the output provides similar spatial resolution to that 
obtained from the more cumbersome, labor-intensive, and 
costly FEM. All local bending and torsional properties of the 
complex structure are then input into computationally effi- 
cient, highly simplified, and accurate analysis methods to 
accurately capture the real physical response of the structure 
without the need for apriori knowledge of the structural 
details of the complex structure. Finally, the computationally 
efficient and highly simplified analysis methods are 
employed during real-time loading of the complex structure 
in order to, obtain real-time, operational loading data. 

This process provides several benefits compared to the 
current methods described above. First, the complexity of the 
FEM can be reduced. And, second, a reduced number of 
degrees -of- freedom are required in this method. These ben- 
efits dramatically reduce the time/money required to accom- 
plish the structural analysis and increase the accuracy of the 
simplified structural models employed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the drawings, which are not necessarily to scale, like or 
corresponding parts are denoted by like or corresponding 
reference numerals. 

FIG. 1 depicts a structure having a plurality of fiber optic 
sensors that comprise a plurality of sections. 

FIG. 2 depicts a cantilever beam with constant cross sec- 
tion under a concentrated load. 

FIG. 3 depicts a cantilever beam with tapered cross section 
subject to tip loading. 
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FIG. 4 depicts a cantilever beam with constant cross sec- 
tion subject to distributed loading. 

FIG. 5 depicts a cantilever beam with tapered cross section 
subject to distributed loading. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

The present invention is an improved method of using 10 
strain data in order to obtain operational loads on complex 
structures. In general, a plurality of fiber optic sensors are 
placed on a complex structure, preferably in a grid type, 
regular pattern. By using fiber optic sensors, a large number 
of sensors may be employed, due to the small size and weight 1 5 
of the sensors. The sensors, in essence, divide the complex 
structure into a plurality of sections. The present invention 
assumes that each of these “sections”, which are defined 
herein as the area from one adjacent fiber optic sensor to 
another, has a consistent structural behavior within each sec- 20 
tion, but the structural behavior may be different in different 
sections. Because the strain is measured at both edges of each 
section, the section lengths can either be selected to be the 
same or may vary. 

After the fiber optic sensors are installed upon the complex 25 
structure, bending and torsional loads can be determined as 
follows. For the bending load, a single, externally applied 
point load, with a known force, is applied to the complex 
structure, normal to the structure’s primary axis. Similarly, 
for torsional loads, a single force couple (moment) at the free 30 
end of the structure, with a known magnitude, is applied to the 
complex, structure. The strain is measured at each fiber optic 
sensor location and, due to the above discretization of the 
complex structure into a plurality of sections, closed-form 
analytical solutions for non-uniform cantilevered beams may 35 
be employed to obtain the bending and torsional structural 
properties of the cross-section of each segment, and, hence, 
the structural properties of the overall complex, structure. 
This part of the method of the present invention obviates the 
need to complete extremely difficult and rigorous finite ele- 40 
ment modeling of the complex structure. 

Next, the complex structure is subjected to operational 
loading. The light weight and small size of the fiber optic 
sensors allows for operational loading without significant 
variance due to sensor equipment. From the structural prop- 45 
erty information for each section and the strain data obtained 
from the sensors during operational loading, the bending 
moment, shear, and loads at each section may be calculated at 
various time intervals selected by the user. 

A more detailed description of the method of the present 50 
invention follows for the bending load. An analogous descrip- 
tion, can also be shown for the torsional load. As an example 
of the method for determining bending load, the basic dis- 
placement equations will be developed for a complex struc- 
ture comprising a uniform cantilever beam with constant 55 
cross sectional properties. The method is built upon the clas- 
sical bending equation for the uniform beam 


d 2 y _ M( x) (1) 60 

TF " ~W 

where y is the vertical displacement, x is the span-wise 
coordinate, M(x) is the bending moment, E is the Young’s 65 
modulus (or elastic modulus), and I is the moment of inertia. 
Equation 1 can be re-written for the non-uniform beam as 


M(X) _ £(*) ( 2 ) 

EI(x ) c(x) 


where €(x) is the strain, c(x) is the distance from a strain 
sensor on the surface of the structure to the neutral axis 
(which is the half-depth for symmetric sections) both 
expressed as functions of the span-wise coordinate, x. 

Referring to FIG. 1, the spatial domain of the complex 
structure 100 is discretized by dividing the structure 100 into 
sections 102 with strain sensing stations i (=0, 1, 2, 3 ... n) 
installed at section junctures x=x z . (e.g. station i=0 at x=x 0 =0 
and station i=n at x=x„ 1 . 

Using a simple structure of length, 1, and half-depth, c(x), 
the sections are defined to be Al(=l/n)di stances apart. 

Equation 2 may be written for each strain, measuring cross- 
section x=x. as 


Mi _ £i ( 3 ) 

Eli ~ ~ 


In which the subscript i indicates the values at the strain 
sensing cross-section x=x z -. 

Or alternatively, if section modulus is defined as 



Cj 


Then equation 3 becomes 


£; 


The first step in the invention is to install a plurality of fiber 
optic strain sensors 104 on the complex structure surface 100 
in such a way as to measure bending strains at desired loca- 
tions. In a preferred embodiment of the invention, the fiber 
optic strain sensors 104 are Fiber Bragg Grating (FBG) sen- 
sors. These sensors are preferred because they are minimally 
obtrusive, ultra-lightweight, easily installed, accurate, 
immune to EMI, and inherently-safe (no joule heating, spark- 
ing). The most preferred configuration and use of FBG sen- 
sors for the present invention employ the Optical Frequency 
Domain Reflectometry (OFDR) technique with hardware 
architecture described in U.S. Pat. Nos. 5,798,521 and 6,566, 
648 which are incorporated herein by reference. This 
approach uses low reflectivity gratings all with the same 
center wavelength and a tunable laser source. The FBGs are 
preferably located on a single optical fiber. This allows hun- 
dreds of strain sensors to be located down the length of the 
fiber; a common configuration is to use 480 FBGs on a single 
20-ft optical fiber spaced at 1 FBG/cm. This configuration 
allows strain measurements to be acquired at much higher 
spatial resolution than other current sensor technologies, 
making it flexible enough to employ a user-selected grating 
density depending on the type of application. 

The number and spacing of the fiber optic sensors 104 
depend upon the overall length of the complex structure 100 
as well as the expected variance of the structural properties of 
the complex structure 100 and may be selected by one skilled 
in the art. As an example, for an aircraft wing, a user may 
employ a fiber optic sensor 104 about every half inch. 
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Referring to FIGS. 1 and 2, the next step is to perform a 
simple test calibration via a single point load of a known 
force. For a beam 200 of length, 1, discretized into x z - sections 
202 each with section lengths of A1 and subjected to a tip load 
of P that places the beam in bending, the flexural rigidity of 5 
the beam can be experimentally determined at each measure- 
ment location, i. Substituting M=P(1-A1), and solving for 
bending stiffness (EI) Z ., equation 3 becomes, 


P(l - iM)Ci (6) 


And subsequently, from equilibrium of forces, the opera- 
tional loads, P z can be determined 

Pi=- AVi ( 11 ) 

The following examples are provided to further illustrate 
the present invention (note that the following examples 
employ finite-element analysis in place of experimental 
results to further demonstrate the invention). Referring to 
io FIG. 2, a cantilever beam with constant cross-section, length 
1=160 in. and subjected to a concentrated tip load, P, of 1000 
lb. is depicted. 


The product of elastic and section modulus in equation 5 
can be determined by calibration at every cross-section by the 1 5 
formula 


{ES)i = 


P(l-iM ) 
£[ 



Table 1 , below, shows the results at each station i obtained 
from a finite element analysis (FEA) and those calculated 
using the method of the present invention. FEA results are 
being used in place of the experimental strains to demonstrate 
the method of the present invention. The table shows excel- 
lent comparison between the flexural rigidity, El calculated 
from the following equation 


In order to obtain operational loads data for the complex 
structure 200, €, is measured while the structure 200 is in 
operation (preferably with the same sensors 204 used during 25 
calibration), and (ES), in equation 6 has been determined via 
a single point calibration, then the bending moment at each 
station, M z , can then be determined for any general set of 
loads applied during operation. Equation 3 becomes 


(EI)i = 


Ebh 3 
12 


(12) 


and from strains generated by FEA input to the equation 6 


30 (£y) = W-Wa (13) 

.. m,s, ( 8 ) ' e ‘ 

Mi = 

c i 


And equation 5 becomes 

M.=(£S).e. (9) 

From the equilibrium of moments, the shear loads. V z , at 
each x z , can then be determined using the following equation 

dM AM/ (10) 

* dx ~ Ax; 


The value obtained from the example analytical “calibra- 
tion” above is then used along with the strains from FEA to 
determine the structure’s moment, shear, and applied loads, 
from equations 8, 10, and 11, respectively. 

Table 1 further shows that the moment from equation 8 
compares well with those computed from M z =P(l-x 1 )=P[l-i 
(Ax)]=P[l-i(Al)]. 

Table 1 also shows that the shear from equation 1 0, and the 
loads from equation 1 1 compares well with those from statics. 


35 


40 


TABLE 1 


r 


Results from simulated experimental results (FEA) and the present invention. 


Calibration 


^ r 


Operation 


A 


index 

distance 

in 

strain 

in/in. 

El calc 
lbs/in2 

El FEA 
lbs/in2 

Moment 

FEA 

Moment 

Calc 

Vi 

FEA 

Vi true 
Calc 

Pi 

FEA 

Pi true 
Calc 

0 

0 

9.59280E-05 

8.33E+09 

8.34E+09 

160000 

160000 

1000 

1000 

0 

0 

1 

20 

8.400 10E-05 

8.33E+09 

8.33E+09 

140000 

140000 

1000 

1000 

0 

0 

2 

40 

7.20000E-05 

8.33E+09 

8.33E+09 

120000 

120000 

1000 

1000 

0 

0 

3 

60 

6.00000E-05 

8.33E+09 

8.33E+09 

100000 

100000 

1000 

1000 

0 

0 

4 

80 

4.80000E-05 

8.33E+09 

8.33E+09 

80000 

80000 

1000 

1000 

0 

0 

5 

100 

3.60000E-05 

8.33E+09 

8.33E+09 

60000 

60000 

1000 

1000 

0 

0 

6 

120 

2.40000E-05 

8.33E+09 

8.33E+09 

40000 

40000 

1000 

1000 

0 

0 

7 

140 

1.19960E-05 

8.33E+09 

8.34E+09 

20000 

20000 

1000 

1000 

0 

0 

8 

160 

6.94000E-07 

8.33E+09 

0.00E+00 

0 

0 

1000 

1000 

-1000 

1000 
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Referring to FIG. 3, a cantilever beam with a tapered cross- 
section 300, length 1=160 in. and subjected to a tip load, P, of 
1000 lb. is depicted. 

Table 2, below, shows the results at each station i obtained 
from the finite element analysis (FEA) and those calculated 5 
using the present invention. The table shows excellent com- 
parison between the flexural rigidity, El calculated from the 
following equation 


(£/); = 


Ebh 3 
12 


10 

(12) 


and from strains generated by FEA input to the equation 6 


TABLE 3 


Results from FEA and the present invention. 


r 


Operation 

A ^ 


Strain FEA 

Moment 

Vi 

Pi 

in./in. 

lb-in 

lb 

lb 

0.00034574 

576671 

7300 

-891 

0.00026909 

448485 

6409 

-809 

0.00020189 

336483 

5600 

-800 

0.00014429 

240483 

4800 

-800 

0.00009629 

160483 

4000 

-800 

0.00005789 

96483 

3200 

-800 

0.00002909 

48483 

2400 

-800 

0.00000989 

16485 

1600 

-776 

0.00000055 

0 

824 

-824 


T\l-iM)d (13) 


The value obtained from the analytical “calibration” above 
is then used along with the strains from FEA to determine the 
structure’s moment, shear, and applied loads, from equations 
8, 10, and 11, respectively. 

Table 2 shows that the moment from equation 8 compares 
well with those computed from M z =P[l-i(Al)]. 


Referring to FIG. 5, a cantilever beam with a tapered cross - 
20 section 500, length 1=160 in. and subjected to a distributed 
load of 800 lb. is depicted. The value obtained from the 
analytical “calibration” in Table 2 is then used along with the 
strains from FEA to determine the structure’ s moment, shear, 
and applied loads, from equations 8, 10, and 1 1 , respectively. 
25 Table 4 shows that the loads from equation 1 1 compares 
well with the applied as shown in FIG. 5. 


TABLE 2 


r 


Results from FEA and the present invention. 


Calibration Operation 

= ,, ^ 


A 


index 

distance 

in 

strain 

in/in. 

El calc 
lbs/in2 

El FEA 
lbs/in2 

Moment 

FEA 

Moment 

Calc 

Vi 

FEA 

Vi true 
Calc 

Pi 

FEA 

Pi true 
Calc 

0 

0 

0.000100269 

8.33E+09 

7.98E+09 

160000 

160000 

1000 

1000 

0 

0 

1 

20 

0.000115283 

5.47E+09 

5.46E+09 

140000 

140000 

1000 

1000 

0 

0 

2 

40 

0.000140685 

3.41E+09 

3.41E+09 

120000 

120000 

1000 

1000 

0 

0 

3 

60 

0.000174987 

2.00E+09 

2.00E+09 

100000 

100000 

1000 

1000 

0 

0 

4 

80 

0.000222254 

1.08E+09 

1.08E+09 

80000 

80000 

1000 

1000 

0 

0 

5 

100 

0.000287898 

5.21E+08 

5.21E+08 

60000 

60000 

1000 

1000 

0 

0 

6 

120 

0.00037428 

2.13E+08 

2.14E+08 

40000 

40000 

1000 

1000 

0 

0 

7 

140 

0.000440201 

6.75E+07 

6.82E+07 

20000 

20000 

1000 

1000 

0 

0 

8 

160 

0.000137267 

1.33E+07 

0.00E+00 

0 

0 

1000 

1000 

-1000 

-1000 


Table 2 also shows that the shear from equation 10, and the 
loads from equation 1 1 compares well with those from statics. 50 TABLE 4 


Results from FEA and the present invention. 


Referring to FIG. 4, a cantilever beam with a constant with 
constant cross-section 400, length 1=1 60 in. and subjected to 55 
a distributed load, P, of 800 lb. is depicted. 


The value obtained from the analytical “calibration” in 
Table 1 is then used along with the strains from FEA to 
determine the structure’s moment, shear, and applied loads, 
from equations 8,10, and 1 1 , respectively. Table 3 shows that 
the loads from equation 1 1 compares well with the applied as 
shown in FIG. 3. 


r 


Operation 

^ \ 


Strain FEA 

Moment 

Vi 

Pi 

in./in. 

lb-in 

lb 

lb 

0.000361084 

576184 

7300 

-899 

0.000369048 

448173 

6401 

-806 

0.000394251 

336284 

5594 

-799 

0.000420633 

240380 

4795 

-800 

0.000445818 

160472 

3995 

-800 

0.000463328 

96561 

3196 

-800 

0.000455266 

48655 

2395 

-801 

0.000369174 

16773 

1594 

-755 

0.000109789 

0 

839 

-839 
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What is described are specific examples of many possible 
variations on the same invention and are not intended in a 
limiting sense. The claimed invention can be practiced using 
other variations not specifically described above. 

What is claimed is: 

1. A method of performing real-time, operational loads 
monitoring of a complex structure, comprising the steps of: 

installing a plurality of fiber optic strain sensors on a sur- 
face of the complex structure; 

applying an external single-point load, of known force to 
the complex structure and obtaining a strain measure- 
ments at each of the fiber optic strain sensors; 

calculating the structural properties of the complex struc- 
ture between each fiber optic strain sensor from the 
strain data obtained at each fiber optic strain sensor; and, 

using the structural properties and the strain measured at 
each fiber optic strain sensor during operational loading 
to obtain real-time operational loads on the complex 
structure. 

2. The method of claim 1, wherein the calculating the 
structural properties step further comprises inputting the 
measured strain of adjacent fiber optic strain sensors into 
closed-form, uniform structure algorithms. 

3. The method of claim 2, wherein the algorithms are 
selected from those to determine bending, torsion, or both for 
the complex structure between adjacent fiber optic sensors. 


10 

4. The method of claim 3, wherein the calculating the 
structural properties step further comprises calculating the 
elastic modulus and a section modulus between each of the 
fiber optic sensors. 

5 5. The method of claim 4, wherein the using the structural 

properties and the strain step further comprises, calculating 
the bending moment between adjacent fiber optic sensors 
during operational loading by multiplying the elastic modu- 
lus by the section modulus from the calculating the structural 
10 properties step and further multiplying the strain from opera- 
tional loading. 

6. The method of claim 5, wherein the using the structural 
properties and the strain step further comprising calculating 
the shear loads between adjacent fiber optic sensors from 

15 equilibrium of the operational loading bending moments. 

7. The method of claim 6, wherein the using the structural 
properties and the strain step further comprises calculating 
the operational loads between the adjacent fiber optic sensors 
from equilibrium of forces of the shear loads. 

8. The method of claim 1, wherein the structural properties 
comprise bending and torsional stiffness. 

9. The method of claim 1, wherein the known force com- 
prises a single-point load for bending. 

10. The method of claim 1, wherein the known force com- 
25 prises an external force-couple for torsion. 



